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Abstract 

Nuclei in dense matter are influenced by the medium. In the cluster mean field approximation, 
an effective Schrodinger equation for the ^4-particle cluster is obtained accounting for the effects of 
the correlated medium such as self-energy, Pauli blocking and Bose enhancement. Similar to the 
single-baryon states (free neutrons and protons), the light elements (2 < A < 4, internal quantum 
state u) are treated as quasiparticles with energies Ea,v{P', T, n n , n p ). These energies depend on the 
center of mass momentum P, as well as temperature T and the total densities n n ,n p of neutrons 
and protons, respectively. No j3 equilibrium is considered so that n n ,n p (or the corresponding 
chemical potentials /J, n ,{j, p ) are fixed independently. 

For the single nucleon quasiparticle energy shift, different approximate expressions such as 
Skyrme or relativistic mean field approaches are well known. Treating the A-particle problem 
in appropriate approximations, results for the cluster quasiparticle shifts are given. Properties of 
dense nuclear matter at moderate temperatures in the subsaturation density region considered here 
are influenced by the composition. This in turn is determined by the cluster quasiparticle energies, 
in particular the formation of clusters at low densities when the temperature decreases, and their 
dissolution due to Pauli blocking as the density increases. Our finite-temperature Green function 
approach covers different limiting cases: The low-density region where the model of nuclear statis- 
tical equilibrium and virial expansions can be applied, and the saturation density region where a 
mean field approach is possible. 

PACS numbers: 21.65.-f,21.45.-v 
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I. INTRODUCTION 



A well established tool in treating many-particle systems is the quasiparticle concept. 
In contrast to free particles, the properties of quasiparticles such as the dispersion relation 
are modified due to the interaction with other particles. This relation between energy and 
momentum is often characterized by an energy shift and an effective mass which depends 
on temperature and density of the medium. A significant part of the interaction can be 
taken into account by introducing this approach. The quasiparticle concept is not restricted 
to elementary particles like nucleons only, but it can also be applied to composed particles, 
i.e. nuclei. (Note that the nucleons itselves are also composed particles.) In this work we 
evaluate these cluster quasiparticle energies within a microscopic approach. 

Nuclear matter is a strongly interacting quantum fluid. In order to treat warm and 
dilute matter (i.e. at subsaturation baryon densities % = n n + n p < 0.16 fm -3 and tem- 
peratures T < 20 MeV) within a systematic quantum statistical approach, we start from a 
nonrelativistic Hamiltonian 

H = y £,E{l)a\a 1 + ± ]T V(12, l'2')a\a\a 2l a v , (1) 

1 12,1'2' 

where {1} denotes momentum HPi, spin cr 1; and isospin T\ characterizing the neutron (n) 
or proton (p) state. The kinetic energy in fermion second quantization a{,a 1 contains 
E(l) = h 2 Pi /2mi, whereas the potential energy contains the matrix element V(12, 1'2') 
of the nucleon-nucleon interaction. 

Since there is no fundamental expression for the nucleon-nucleon interaction (like, e.g., 
the Coulomb interaction in charged particle systems [l!]), it is taken to reproduce empirical 
data such as the nucleon scattering phase shifts. Different parametrization are in use. Simple 
potentials as proposed by Yukawa, Yamaguchi, Mongan, Gogny and others are based on two- 
nucleon phase shifts and can be used for exploratory calculations. For detailed calculations 
one can use more sophisticated potentials such as PARIS and BONN or their separable 
representations [2| • To obtain the empirical parameter values of nuclear matter at saturation 



density, three-body forces have 



Argonne AV18/UIX potential [3] has been used to calculate light nuclei 

Near saturation density, nucleons can be treated as quasiparticles. Semiempirical ap- 
proaches such as the Skyrme contact pseudopotential [5f] and relativistic mean field (RMF) 



?een introduced in the Hamiltonian ([T]). In particular, the 



approaches [6j parametrize the quasiparticle shift as a function of densities and tempera- 
ture. Microscopic approaches such as Dirac-Brueckner Hartree Fock (DBHF) 7] give an 
appropriate description of the thermodynamic properties of warm and dense matter. For a 
recent review of the nuclear matter equation of state (EoS) see Ref. jg]. 

In the low-density limit, nuclear matter at finite temperature is a mixture of free nucle- 
ons and nuclei in chemical equilibrium as described by a mass action law. This chemical 
picture, where bound states are treated as new species, is also denoted as nuclear statistical 
equilibrium (NSE) and should be recovered as the low-density limit of a quantum statistical 
approach to nuclear matter at finite temperatures. 

At increasing densities, the simple NSE approach becomes invalid because of interactions 
between the nucleons and nuclei. A systematic coherent description should include scattering 
states, avoid double counting of many-particle effects, and avoid semiempirical concepts such 
as excluded volume. This can be achieved within a physical picture where bound states 
are produced by the interaction, investigating the few-particle propagator in the medium. 
The chemical picture can be used as a guideline to select important contributions within a 
perturbative approach and performing partial summations of the corresponding Feynman 
diagrams. In the case of charged particle systems, where the interaction is given by the 
Coulomb potential, the use of the chemical picture to introduce bound states in the partially 
ionized plasma has been extensively worked out, see Ref. 

The in-medium wave equation for a system of A nucleons is derived in Sec. HI1 Bound 
states describe nuclei with energy eigenvalues EA,v(P;T,n p ,n n ) where A denotes the mass 
number. The index v specifies the internal state of the A-nucleon system such as spin, isospin 
and excitation. In a homogeneous system, the center-of-mass momentum P of the cluster 
is conserved and can be used as quantum number. Assuming (local) thermal equilibrium, 
due to the influence of the surrounding matter these energy eigenvalues will depend on 
three parameter: The temperature T and the total number densities n p , n n of protons and 
neutrons, respectively [or the baryon density % = n n + n p and the asymmetry parameter 
a = (n n — n p )/nB}- We do not consider /3-equilibrium due to weak interaction processes. 
Explicit results for the corresponding quasiparticle shifts are given for 2 H (deuteron d), 3 H 
(triton t), 3 He (helion h), and 4 He (a-particle) in Sec. IIVI As an application, the generalized 
Beth-Uhlenbeck formula is outlined in Sec. |V] and the nuclear matter EoS is considered 
which contains the quasiparticle energies of the single-nucleon states as well as of the light 
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nuclei. 



II. MANY PARTICLE APPROACH 

A. Single-particle spectral function and quasiparticles 

The basic equations of a Green function approach to the many-nucleon system can be 
found in different textbooks and papers, see js, Here we give only some final results. 

The propagation of a single- nucleon excitation (a\(t)a lt ) can be expressed in terms of the 
single-particle spectral function Si(l,u) as 

(a\(t)a v ) = 5 1>v J ^e**fa(cj)Si(l,u>) , (2) 

where 

f A)Z (u) = [exp(/3(o; - -{A- Z)fi n )) - (-1) V (3) 

is the Fermi or Bose distribution function which depends on the inverse temperature (3 = 
1/(/cbT) and the chemical potentials \i v ,\i n (instead of the isospin quantum number T\ we 
use the charge number Z). A special case is the single-particle density matrix (a\a v ) which 
can be used to evaluate the equation of state (EoS) for the nucleon density 

d 3 P 

n T (/?,/i p ,/i n ) = — \ ](a[a 1 )5 TjTl = 2 / — — - / ^fi jZ {u)Si{\, u) , (4) 



where Q is the system volume, and summation over spin direction is collected in the factor 2. 
Both the Fermi distribution function and the spectral function depend on the temperature 
and the chemical potentials n P , fi n not given explicitly. We work with a grand canonical 
ensemble and have to use the EoS (j4j) to replace the chemical potentials by the densities 
n p , n n . For this equation of state, expressions such as the Beth-Uhlenbeck formula are 
discussed below in Sec. IVl 

The spectral function is related to the self-energy according to 

9n ) = 2ImS(l,o;-»0) 

U ' ' \uj-E{\) -ReE(l,a;)] 2 + [ImS(l, - «0)] 2 ' U 

where the imaginary part has to be taken for a small negative imaginary part in the fre- 
quency. The solution of the relation 

£?! qu (l) = E(l) + ReS[l, #ni)] (6) 




T2 : ex 








T^ex 




+ > 




+ 





+ 



FIG. 1: Cluster decomposition of the single- nucleoli self-energy. The index 'ex' denotes full anti- 
symmetrization including all exchange diagrams. For bound states with A < 4 the direct ladder-T 
matrices can be taken, without the full exchange at the end, as long as the particles have different 
internal quantum numbers a, r. The Fock term results from T2 )ex (first order in V, exchange term, 
see also App. |B|) , 

defines the single-nucleon quasiparticle energies -E^ u (l) = E(l) + A£' SE (1). Expanding for 
small Im£(l,z), the spectral function yields a 5-like contribution, and the densities are 
calculated from Fermi distributions with the quasiparticle energies so that 

< u (/3,^,^) = ^/i,z[^r(i)] 



(7) 



Pi 



follows for the EoS in mean field approximation. This result does not contain the contri- 
bution of bound states and therefore fails to be correct in the low-temperature, low-density 
limit where the NSE describes the nuclear matter EoS. 

As shown in Refs. [lo|, llj], the bound state contributions are obtained from the poles 
of ImS(l, z) which cannot be neglected expanding the spectral function. A cluster decom- 
position of the self-energy has been proposed, see Fig. [TJ The diagrams are calculated 
as 

£(M,) = E ]T Gf_ 1 (2,...,A,n x -z u )T A (l...A,l'...A',n x ). (8) 

A>1 ft^,2...A 

The free (A — I) quasiparticle propagator G^_ 1 and the T^ matrix are given in App. \K\ The 
Ta matrices are related to the A-particle Green functions which read in bilinear expansion 

1 



G A (1...A, I'... A', z A ) = £ V^p(l ...A) 



AvP 



A' ...A'). 



(9) 



The A-particle wave function tpAvp(l ■ ■ ■ A) and the corresponding eigenvalues E A \(P) result 
from solving the in-medium Schrodinger equation, see the following Sec. Ill BL Besides the 
bound states, the summation over the internal quantum states v includes also the scattering 
states. 
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The evaluation of the equation of state in the low-density limit is straightforward, see 
App. [A] Considering only the bound-state contributions, we obtain the result 



1 . 

(T,// p ,/v) = — ^ Z f AiZ [Ef iU (P;T, fi p , fi n )] , 

A,u,P 

(T )/Wn ) = i J2( A ~ Z )fA,z[ET, u (P;T,fi p ,fi n )} (10) 



A,u,P 

for the EoS describing a mixture of components (cluster quasiparticles) obeying Fermi or 
Bose statistics. The NSE is obtained in the low-density limit if the in-medium energies 
E^ U (P;T, fi p , fi n ) can be replaced by the binding energies of the isolated nuclei E Au (P). 
Note that at low temperatures Bose-Einstein condensation may occur. 

We discuss the cluster decomposition of the single-nucleon self-energy, Fig. [TJ in com- 
parison with other approximations. Restricting the cluster decomposition only to the con- 
tribution of two-particle correlations, we obtain the so-called ToG approximation. In this 

nn 

approximation, the Beth-Uhlenbeck formula is obtained for the EoS, as shown in [10, 
The relation between this T matrix approach and the Brueckner G matrix approach was 
discussed in detail in Ref. 12j. Extended work has been performed using sophisticated 
interaction potentials to evaluate the quasiparticle energies in the DBHF approximation, for 
recent reviews see Refs. 0, H, fl4[. 

Replacing the T2 matrix in Born approximation with the interaction potential V, we 
obtain the Hartree-Fock approximation 

A HF (1) = J>(12, 12) - V(12, 21)]/ liT2 [£(2) + A HF (2)] . (11) 

2 

In this approximation, all correlations in the medium are neglected. The self-energy does 
not depend on frequency, i.e. it is instantaneous in time, with vanishing imaginary part. 



B. Effective wave equation for the yl-nucleon cluster 

We consider the propagation of an A-nucleon cluster in warm, dense matter which is 
described by the A-particle Green function Ga- The solution will not only enter the cluster 
decomposition of the single-nucleon self-energy as considered in Sec. Ill Al to calculate the 
contribution of bound states to the nuclear matter EoS in a systematic way. It determines 
also the cluster decomposition of other quantities such as the polarization function, the 
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dynamical structure factor etc. We can proceed as above in the single-nucleon case and 
investigate higher order correlation functions. The A-nucleon spectral function Sa can be 
introduced which is related to G a- Cluster quasiparticle excitations are determined by 5-like 
peaks in the spectral function Sa- 

For the A-particle Green function, the perturbation expansion can be represented by 
Feynman diagrams. New elements such as the A-particle self-energy can be introduced, see 



Ref. 



151 ] . In the low-density limit, the bound and scattering states of the A-particle 



cluster are obtained performing the partial summation of ladder diagrams. 

The interaction of the A-nucleon cluster with the surrounding nucleons can be considered 
in the cluster-mean field approximation given in App. [B] Since nuclear matter can form 
clusters, the interaction of the A-nucleon cluster under consideration is taken with an arbi- 
trary cluster B of the surrounding matter in first order of the nucleon-nucleon interaction, 
but full antisymmetrization of all nucleons in the cluster A and B. As a result, not only 
the Hartree shift due to the interaction with the surrounding nucleons in free as well as in 
cluster states is obtained. Also the the Pauli blocking terms due to the occupation of the 
phase space by free nucleons as well as nucleons which are bound in clusters is taken into 
account. In higher orders, the Bose-enhancement or Pauli blocking of nucleonic clusters is 
obtained. 

At present, the full self-consistent solution of the cluster- mean field approximation is out 
of reach. Only special cases can be solved such as a cluster nuclei where an a cluster is 
considered in a surroundings consisting of a few a clusters, with full antisymmetrization on 
the nucleonic level. 

We calculate the modification of the /L-nucleon cluster due to the surroundings consid- 
ering only self-energy and Pauli-blocking terms in the cluster-mean field approximation. 
Correlations in the medium are neglected. In momentum space, the distribution of free 
nucleons, forming a Fermi sphere, and the bound states which are characterized by a wave 
function, are different functions, but the total amount of the occupied phase space is deter- 
mined by the total density of the protons or neutrons. Considering the uncorrelated medium, 
the phase space occupation is described by a Fermi distribution function normalized to the 
total density of nucleons, see App. [HI 

1(1) = - « ZZ. (z^L. e -ET(l)/T (12) 

Jl{ ) exp[E* u (l)/T-fr T /T) + l 2 \mT ) { ' 
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in the low-density, non-degenerate limit (/2 T < 0). The effective chemical potential p, T is 
determined by the normalization condition 2Q~ 1 Y2 p fi{p) — n r, where r denotes isospin 
(neutron or proton), and has to be expressed in terms of these densities. 

In ladder approximation, the A-particle Green function obeys a Bethe-Salpeter equation 
(BSE) 

G A (1...A, I'... A', z A ) = Gf (1...A z A )5 lv . . . 5 AA , 

+ Gf(l...A,z A )V A (l...A,l"...A")G A (l"...A",l'...A',z A ), (13) 

1"...A" 

where V A (1...A, 1' . . . A') = J2i<j V^Oi i'j') Ylk^ij ^ kk ' ls ^ ne interaction within the A-particle 
cluster. z A denotes a fermionic or bosonic Matsubara frequency. The free A-quasiparticle 
Green function results as 

V- A > ZA) - z A -ET{l) E?(A) ' [ 1 

The solution of the BSE is given by the bi-linear expansion (Q. The A-particle wave 
function and the corresponding eigenvalues follow from solving the in-medium Schrodinger 
equation 

[ET(1) + " " " + ^) - EZ(PMa,p(1 ...A) 

+ ^[l-fi(*)-hUW(ij,i'f) H5 kk ^ AvP (l'...A') = 0. (15) 

1'.. .A' i<j k^i,j 

This equation contains the effects of the medium in the quasiparticle shift as well as in the 
Pauli blocking terms. Obviously the bound state wave functions and energy eigenvalues 
as well as the scattering phase shifts become dependent on temperature and density. Two 
effects have to be considered: the single-nucleon quasiparticle energy shift and the Pauli 
blocking. 

Detailed investigations of Eq. (j!5p have been performed for A = 2, see Ref. |10] . de- 
scribing interesting physics. From the solution of this in-medium two-particle Schrodinger 
equation or the corresponding T matrix, the medium dependent scattering and possibly 
bound states are obtained, as well as the formation of a quantum condensate, including the 
crossover from Bose-Einstein condensation to Cooper pairing. Due to the self-energy shifts 
and the Pauli blocking, the energy of the deuteron, E^ n (P; T, fi p , /i n ), as well as the scattering 
phase shifts, S T (E, P;T, fi n ) (r denoting the isospin singlet or triplet channel), will de- 
pend on the temperature and the chemical potentials. For a separable interaction V(12, 1'2') 
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like the PEST4 potential [2j, solutions of the in- medium two-particle Schrodinger equation 



can be found, e.g., in Refs. [1 
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17| . We will evaluate the medium shift of the binding 



energy in the following section using perturbation theory. 



III. EXPRESSIONS FOR THE LIGHT NUCLEI QUASIPARTICLE SHIFTS 
A. Solution of the few-nucleon effective wave equation 

The quasiparticle energies in warm and dense nuclear matter, E A \ t (P;T,n p ,n n ), are well 
defined functions of temperature and total proton and neutron densities, given by a peak in 
the v4-particle spectral function. To evaluate these quantities for infinite matter from a first- 
principle quantum statistical approach, we have to perform some approximations. Note 
that the quasiparticle shifts can also be introduced phenomenologically fitting empirical 
data. This is well-known for the single-nucleon quasiparticle energy which can be adapted 
to reproduce the structure of nuclei p]. 

The solution of the few-body in-medium Schrodinger equation fll5D for separable interac- 



tion is simple in the case A=2. For A=3, 4, a Faddeev approach can be used, see [18|, [19]. 
To obtain explicit expressions for the quasiparticle energy shifts, we will apply perturbation 
theory, which can be justified in the low-density region. Denoting the unperturbed wave 
function of the A-nucleon cluster with v?Aiap(1 ■ ■ ■ A), we have 

EZ(P) ~ E • • • • • • A 1' • • ■ A') VAvP {l' ...A'), (16) 

l...A,V ...A' 

where the form of H eS (l . . . A, V . . . A') is given by Eq. ([15!) after symmetrization. 

Before calculating the in-medium quasiparticle energy eigenvalues ffl6|) . A < 4, we 
first consider the isolated A-nucleon problem to determine the unperturbed wave function 
ifAvpi)- ■ ■ ■ A). Extended analysis has been carried out by Wiringa et al. using Green's func- 
tion Monte Carlo calculations 4|, where the AV18/UIX potential was used. Some properties 
are given in Tab. [B There is excellent agreement between theory and experimental data 

In contrast to the second virial coefficient of the EoS given below in Sec. EJ which is 
determined by on-shell properties (binding energy and scattering phase shifts of the two- 
nucleon problem) determined from experiments, the nucleon interaction potential and the 
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TABLE I: Light cluster properties at zero density 





binding 


mass 


spin 


rms-radius 


rms-radius 




energy 






(charge) 


(point) 




[MeVl 

L J 


[MeV/c 2 l 

L /J 




[fm] 

L J 


[fm] 

L 1 


n 





939.565 


1/2 

/ 


0.34 





p(!H) 





938.783 


1/2 


0.87 





d ( 2 H) 


-2.225 


1876.12 


1 


2.14 


1.96 


t ( 3 H) 


-8.482 


2809.43 


1/2 


1.77 


1.59 


/i ( 3 He) 


-7.718 


2809.41 


1/2 


1.97 


1.76 


a ( 4 He) 


-28.30 


3728.40 





1.68 


1.45 



bound state wave function must be known to evaluate the quasiparticle energies . This 
ootential is derived fitting empirical data which have to be reproduced. In previous work 
lo| . the PARIS potential [2j has been used. Alternatively, the AV18/UIX potential 4] 
can be taken. Here, we will use simple separable interaction potentials, which are fitted 
to the binding energy and the rms radius value of the respective cluster, in order to obtain 
analytical expressions for the quasiparticle shifts as function of P, T, n p , n n . Details are given 
in App. O 



1. Gaussian wave function approach 

For the separable pair interaction in the A-nucleon cluster a Gaussian form is taken, 

V V {12, 1'2') = \ v e e ^ , (17) 

where v = {A, Z} := {d, t, h, a} denotes the cluster under consideration. 

To solve the A-particle Schrodinger equation, a variational approach is used. Two dif- 
ferent classes of functions are considered. First, a Gaussian wave function is taken which 
allows for analytical expressions in evaluating the shifts. In the following Sec. IIII A 21 a 
better, but more complex Jastrow ansatz is considered which, however, allows in general 
only for numerical evaluation of the shifts. 

For Gaussian wave functions, the center of mass motion can be easily separated. For 
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TABLE II: Light cluster wave function parameter at zero density from the Gaussian approach 



V 


K 








[MeV fm 3 ] 


[fm" 1 ] 


[fm- 1 ] 


d ( 2 H) 


-3677.2 


0.753 


0.625 


t ( 3 H) 


-1670.0 


1.083 


0.889 


h ( 3 He) 


-1957.5 


0.960 


0.804 


a ( 4 He) 


-1449.6 


1.152 


1.034 



vanishing center-of-mass motion, P = 0, we have 

7 .-:;""-(/',.../m! = — !— e-^ + -^)/Bi5 (18) 
norm,, 

with the normalization Y2 Pl PA \<Pv nss {Pi- ■ -Pa)\ 2 — 1- The parameter B v is fixed by the 
nucleonic point rms radius a/ (r 2 ) u , given in Tab. [H 

ff 2_ 3(A-l) 

see App. Values for are presented in Tab. ILT1 

Next we are interested in the parameter values X u , % of the potential (JTTj) which yield the 
binding energy of the nucleus as well as the nucleonic point rms radius, given in Tab. 
[H using the variational ansatz (fl~8l) . Calculating the kinetic energy as well as the potential 
energy for a Gaussian wave function with range parameter B u , we have 
(0) _ 3(A-1) A(A-1) !lBl 

E » -^^ + ^^kTW' ( } 

Varying B u , the potential parameters which have the minimum energy at B v , consistent 
with the rms radius value, are shown in Tab. (TTl 



2. Jastrow wave function approach 

Of course, the variational solution of the separable Gauss potential (TTT|) by using Gaus- 
sians is not optimal, which is clearly seen for the two-nucleon system where the exact solution 
is known. A Jastrow ansatz which reproduces this exact solution for A = 2 is given by 



1 p 4a^ 

*>«~(fi...ao = -ji—^-. (2i) 



i<3 462 
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TABLE III: Light cluster wave function parameter at zero density from the Jastrow approach 



V 


K 




a u 


o v 




[MeV fm 3 ] 


[fm" 1 ] 


[fm- 1 ] 


[fm" 1 ] 


d ( 2 H) 


-1287.4 


1.474 


1.474 


0.2317 


t ( 3 H) 


-1283.8 


1.259 


1.742 


0.592 


h ( 3 He) 


-1527.0 


1.105 


1.528 


0.542 


a ( 4 He) 


-1272.9 


1.231 


2.151 


0.912 



The prefactor N v is determined by the normalization condition. The results for the pa- 
rameter values X u , j u of the interaction potential, which reproduce the binding energies and 
rms values of the cluster (see Tab. I), as well as the parameter values a u ,b u characteriz- 
ing the wave function (1211) within the variational approach, are given in Tab. IIHI For the 
calculations see App. O 



B. Expressions for the quasiparticle shifts 



1. Nucleon quasiparticles 

The in- medium single- nucleon dispersion relation Ef u (P) can be expanded for small mo- 
menta P as 

E ^(P) = —P 2 + AEf E (P) = AE 1 SE (0) + -^P 2 + 0(P 4 ) , (22) 

where the quasiparticle energies are shifted by AEf E (0), and m\ denotes the effective mass 
of neutrons {t\ = n) or protons (n = p). Both quantities, APf E (0) and m*, are functions 
of T,n p ,n n characterizing the surrounding matter. 

Different expressions are used to parametrize the nucleon quasiparticle shift at subsatura- 
tion density. The Skyrme parametrization 3] of AEf E (0), see App. ICl is used in a standard 



approach to the nuclear matter EoS by Lattimer and Swesty [21[. Different improvements 
have been performed to optimise the calculation for various nuclei. Alternatively, relativistic 
mean-field (RMF) approaches have been developed starting from a model Lagrangian which 
couples the nucleons to mesons. The relativistic quasiparticle energy is given by 



E? U (P) = yj [m T c 2 - S(n n , n p , T)] z + K^P 2 + V T (n n , n p , T) . (23) 
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Expressions for S and V T can be given 22j determining the nucleon quasiparticle shift and 
the effective mass, which aim to obtain the nuclear matter EoS for supernovae collapses. 
Recent work on RMF parametrization intend to reproduce properties of nuclei, but 

are also in agreement with microscopic DBHF calculations. 

The EoS of asymmetric nuclear matter has been investigated in the low-density region 
below the nuclear saturation density [ljj], and expansions of the quasiparticle shift in powers 



of the density have been given. The Dirac mass and the Landau mass are considered [13]. 
A more detailed discussion of the nucleon quasiparticle approach which is very successful in 
describing nuclear matter near saturation density cannot be given here. 



2. Deuteron quasiparticles 

In the low-density limit, the shift of the deuteron binding energy can be calculated from 
the in-medium Schrodinger equation ( fl5|) taking the medium modifications as correction. 
Within perturbation theory, the shifts of the solution of this medium modified wave equation 
are given by the self-energy term and the Pauli blocking term, 

ET(P) = ET + ^-P 2 + A£f(P) + A£j auli (P) , (24) 

ma ~ 2m is the deuteron rest mass. After separation of the center-of mass motion, the 
free deuteron wave function <p d {qi) is taken as a function of the relative momentum q\ = 
(p2 —pi)/2 (neglecting the proton-neutron mass difference). Perturbation theory gives for 
the self-energy term due to the single- nucleon self-energy shift AEf E (P) = El u (P) — E{1) 
in the medium part of the effective Hamiltonian (Tl~5"j) 

A£f(P) = ^5>^i) 

Nd = X]<ji IVdiqi)^ is the normalization. For a separable interaction, the Pauli blocking 
term reads 

A £Pauli (p) = _ J_ J2 V *M) 
<?i,9i' 

where the distribution function fi(p) is given by Eq. (fl2"|) . After applying the unperturbed 
Schrodinger equation, the interaction in Eq. (1261) can be eliminated. With the reduced mass 
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(25) 



fn 



P 



+ qi)+ f P 



V(qi,qi')<p d {qi), (26) 



m/2 we obtain 



(27) 



We give some more explicit expressions. Similar to the effective mass representation for 

the single-nucleon quasiparticle states, we can introduce the deuteron shift and effective 
mass according to 

ET(p) = E (o) + + _V_ p2 + 0(p4) (2g) 

The self-energy contribution to the deuteron shift contains the contribution of the single- 
nucleon shift 

AE rigid shift (Q) = A £SE (Q) + A £f(0) . (29) 

The rigid shift of the nucleons is also present in the continuum states and will not change 
the binding energy of the bound states. It can be incorporated in the chemical potentials 
/i n ,/x p , similar to the rest mass of the nucleons. 

If we consider the effective mass of the nucleons m* ^ m T , we get the contribution 

ABf-W = ±- E | w(a)P * + ^ - 2) (rf + £) . ,30, 

The first term contributes to the deuteron shift AP| E (0), the last term (oc P 2 ) to the 
deuteron effective mass. With Am^ = m* d — we find for the self-energy contribution 



Amf 1 



m p m r , 



R + ■ (31) 

For the further evaluation, we need the free deuteron wave function. For the solution (pTj) 
of the Gaussian interaction ( [171) . used in the Jastrow ansatz with parameter values given 
in Tab. IHIl we obtain for the Pauli blocking shift ( 1271) the expression APj aull,Jastrow (P) = 
(n p + n n )£Pj aull ' J (P) with 

1 {2tvH 2 \ 3/2 f , „ e" 2g2 / a ' 



^(P).i^ y d „ 2 _i_i_ e -*^ + o (32) 



x — e 2mT — e 2mT — qi — E, ■ - - 

k 2 Pq V / W ' X 



^(-l + e^^(l + 4|)erfc[^] 

a d 
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taking the nucleons in the medium as non-degenerate what, at finite temperatures, is justified 
in the low-density limit. (In the degenerate case, the Boltzmann distribution has to be 
replaced by the Fermi distribution.) In particular, we find the quasiparticle shift at P = 



5E 



Pauli,J 



mT 



3/2 tfal 



2m 



1 I n a d 



-1/2 



erfc 



4mT 



(33) 



erfc 



V2 



Values for the deuteron quasiparticle shift are given below in Tab. IIVI 

To obtain the Pauli blocking contribution to the deuteron effective mass, we can expand 
the Pauli blocking shift ( 1321) for small values of P. We will not give the corresponding 
expressions here. The effective mass approximation for the deuterons in warm dense matter 
is limited to small values of P. For arbitrary P, after averaging over the direction between 
(fi and P, we get the approximation 



AE 



Pauli, Jastrow 



(P) » AE, 



Pauli, Jastrow 
d 



(0)e- 



(34) 



Together with Eqs. (1531 . (I30l . (|29|) . we obtain the momentum-dependent quasiparticle 
energy ( j24j) describing the deuteron in matter. 

In contrast to the self-energy shift, the Pauli blocking shift depends strongly on tempera- 
ture by the following reasons. The binding energy per nucleon characterizes the extension of 
the wave function in momentum space. The form of the Fermi distribution function is deter- 
mined by the temperature T which we consider here to be of the same order as the binding 
energy per nucleon. Further, the Pauli blocking term is strongly depending on the center 
of mass momentum P because of the overlap of the deuteron wave function in momentum 
space with the Fermi sphere. Therefore, the bound states with high momentum P are less 
modified by the Pauli blocking effects. The evaluation of the deuteron quasiparticle energy 
is further improved using a better wave function, an appropriate interaction potential and 
avoiding perturbation expansions. Comparison with other potentials such as the Yamaguchi 
Lorentzian formfactor has been performed which give only small changes (below 10 %) to 
the deuteron quasiparticle shift. 
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3. Tritium and Helium in matter 



We now consider the clusters with A=3 (t, h) and A=4 (a). Some of the relations given 
for the deuteron case in Sec. IIIIB2I can be generalized to higher values of A. In the low- 
density limit, the shift of the cluster binding energy can be calculated from the effective 
Schrodinger equation (fl5l) taking the medium modifications as correction, 

E?(P) = + ^—P 2 + AE„(P) (35) 

with v = {A, Z} = {t, h, a}, and m v « Am denoting the rest mass of the cluster. Within 
perturbation theory, the shift AE I/ (P) = AE^ E (P) + AP^ aull (P) consists of the self-energy 
and Pauli-blocking term which can be calculated from the effective wave equation (fl5l) with 
the unperturbed wave function ip u (pi, . . . , Pa) • After separation of the center-of mass motion, 
(p u (qi, ■ ■ ■ , Qa-i) is a function of the remaining Jacobian momenta qi, see App. O We find 
for the energy shifts due to the single-nucleon self-energy shift in the wave equation ( fl5j) 



A£f (P) = ^ E l^)! 2 i AE n E (Pi) + ■■■ + ^(p A )} . (36) 
and for the Pauli blocking term 

AE Pauii (p) = J- V *M) [fi(pd + VipiViM) II 5 vu,v' h fMi) , (37) 

V qi,qi' i<j k^i,j 

where N„ = Ylq \^v{Qi)\ 2 ls the normalization. The single-nucleon momenta pi have to be 
expressed by Jacobi momenta <fj including the center of mass momentum P. The Pauli 
blocking term results as where also the momenta pi have to be expressed in terms of the 
Jacobi momenta q}. 

As for the single-nucleon and deuteron states, we can introduce the quasiparticle shift 
and effective mass according to 

ET{P) = Ef + A£„(0) + ^P 2 + 0{P A ) . (38) 

The self-energy contribution to the quasiparticle shift contains the contribution of the single- 
nucleon shift 

frigid shifty = ^ _ Z ) A £SE( ) + ZA £SE(0) . ( 39 ) 
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Evaluating the distribution functions fA,z[E^ v (P)], Eq. (fTUl) . this contribution to the quasi- 
particle shift can be included renormalizing the chemical potentials \x n ,\i v . The further dis- 
cussion of the self-energy contribution to the quasiparticle shift and effective mass can be 
performed in analogy to the deuteron case and will not be repeated here. 

Whereas in the deuteron case it was possible to eliminate the interaction potential calcu- 
lating the Pauli blocking shift, Eq. (1271) . this is not longer possible in the case A = 3, 4, and 
we have to evaluate expression fl37j) . Analytical expressions for the Pauli blocking shift will 
be given in the following Sec. IIV} based on the more simple Gaussian ansatz for the unper- 
turbed wave function of the cluster. For the Jastrow ansatz (1211 . results for ( l >E' paull ' J (0) are 
given below in Tab. IIVI 

As already discussed for the deuteron case, after angular averaging the Pauli blocking 
shift can be approximated as 

2 2 

AP Pauli (P) « AP Pauli (0) e~^r . (40) 

It leads to the dissolution of the bound states below saturation density, starting at the 
Mott baryon density n^ Mott (T, a) depending on temperature T and asymmetry parameter 
a, where the A-nucleon bound state with P = merges in the continuum of scattering states 



IV. RESULTS FOR THE CLUSTER QUASIPARTICLE SHIFTS AT LOW DEN- 
SITIES 

We collect some results for the quasiparticle energy 

E?(P; T, n n , n p ) = Pf + + APf (P) + A^ Pauli (P) + AE^\P) (41) 

for the light elements v = {d, t, h, a} = { 2 H, 3 H, 3 He, 4 He} in warm, dense nuclear mat- 
ter. In addition to the single-particle self-energy shift AE$ E (P) and the Pauli blocking 
term AP^^^P), the quasiparticle energy shift contains also the Coulomb shift AP^ oul (P) 
which will not be elaborated here. The remaining two contributions are calculated within 
perturbation theory. 

The quasiparticle self-energy shift AP^ E (P) is caused by the self-energy shift APf E (P) 
of the single-nucleon energies. For the A-nucleon system, it contributes to the bound state 



17 



energies (nuclei) as well as to the energy of scattering states, in particular the edge of the 
continuum of scattering states, see Ref. In addition to the rigid shift, the quasiparticle 
self-energy shift at zero momentum A^ E (0) = A££ igidshlft + 3(A - l)h 2 Bl(m - m*)/(8m 2 ) 
contains the contribution of the effective nucleon masses, calculated for Gaussian wave func- 
tions where B u is given by Eq. (j!9p . Similarly, the self-energy contribution to the cluster 
effective mass can be calculated. 

The Pauli blocking contribution AE^ auh (P) may be taken in the approximation fHUj) 
so that we discuss only the cluster quasiparticle shift AE^ aull (0) here. The Hamiltonian 
describing the unperturbed cluster contains parameters for the interaction potential which 
are determined such that the correct binding energy and point rms radius of the nuclei, 
given in Tab. [U are reproduced. Two different approximations are considered. 

(i) The variational function to solve the isolated cluster Schrodinger equation is taken as 
a Gaussian wave function with the parameter B v . A Gaussian potential (fl7|) is constructed 
which reproduces the correct values for the binding energies and the cluster point rms radii. 
The corresponding values for \ v ^ v are given in Tab. [Ill Now, having the interaction and 
the wave function to our disposal, we can calculate the shift of the binding energies in the 
low-density region using perturbation theory. Expanding with respect to the densities n n , n p , 
we have 

AjE Pauli (0) = [(1 + + (1 _ ^) np ] 5 £PauU (T) + (43) 

with ad = 0, a t = 1/3, ah = —l/3,a a = 0. In the non-degenerate case (/i r < 0), expressions 
for 5E^ auh (T) are obtained in analytic form as 

^Pa U n,G (T) = A(A - 1) (2*h* V /2 KBlt 



47T3/2 V mT ) { Bl + 273)3/2 (jB 2 + 2^2 + g|| (jB 2 + d ^2 )f ,2 ' 

(43) 

where q = 2, c t = Ch = 24, c Q = 16; dd = 0, d t = d h = 14, d a = 10. 

(ii) The same procedure as in (i), only the Jastrow function (j2~T|) is taken as variational 
ansatz for the wave function. The Gaussian potential (TT7|) with parameter values Xu,1u, see 
Tab. Hm reproduces the correct values for the binding energies and the cluster point rms 
radii. This variational ansatz gives the exact solution (I3"3"|) of the two-nucleon problem, so 
that it is expected to yield better results also for the higher clusters. However, no analytic 
expressions are obtained for A > 2. 
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TABLE IV: Temperature dependence of the first order Pauli blocking shift for A = 2, 3, 4. Com- 
parison between the Gauss (G) and the Jastrow (J) approach. T in [MeV], 5E^ auh in [MeV fm 3 ]. 



T 


^Pauli,G 


^Pauli.J 

a 


^Pauli.G 


5 pPauli,J 


^Pauli.G 
" h 


^Pauli.J 
' h 


^Pauli,G 


^Pauli.J 


20 


157.9 


172.9 


482.7 


470.77 


438.9 


432.52 


967.7 


950.16 


15 


228.1 


235.4 


651.7 


628.95 


604.0 


588.70 


1263.8 


1237.72 


10 


371.4 


352.5 


967.7 


905.37 


906.7 


871.81 


1749.1 


1711.83 


9 


418.8 


389.2 


1037.9 


986.54 


998.4 


957.23 


1884.2 


1844.90 


8 


477.1 


433.3 


1140.0 


1081.19 


1106.6 


1058.08 


2037.5 


1996.84 


7 


550.3 


487.7 


1260.0 


1192.82 


1236.0 


1178.72 


2212.9 


2171.86 


6 


644.4 


556.2 


1402.7 


1326.29 


1392.6 


1325.31 


2415.0 


2375.56 


5 


768.8 


645.2 


1574.4 


1488.50 


1585.4 


1506.83 


2649.8 


2615.75 


4 


939.0 


766.2 


1784.3 


1689.66 


1827.2 


1736.97 


2925.4 


2903.94 


3 


1182.8 


941.0 


2045.4 


1945.59 


2137.3 


2037.73 


3252.1 


3259.11 


2 


1553.9 


1212.0 


2377.0 


2282.48 


2546.2 


2447.21 


3644.7 


3720.67 


1 


2169.8 


1756.9 


2808.9 


2747.92 


3104.5 


3038.96 


4123.0 


4434.30 



In Tab. [IV] the comparison between the Gaussian (8E„ ulh ) and the Jastrow ansatz 
(&E* >aull ' J ) is shown for different temperatures. The differences are small (below 5 %) for 
A = 3,4 so that the analytical expression (j4"3"l) can be used. The differences are larger for 
A = 2 because the deuteron wave function is not well approximated by a Gaussian. However, 
for A = 2 we can take the Jastrow result which is given analytically in Eq. ( J33l) . 

The wave functions of the isolated clusters determine the quasiparticle shifts within per- 
turbation theory, in particular the Pauli blocking shift. The correct reproduction of the 
rms radii of the nuclei is important in order to estimate the region in phase space, which is 
needed to form the bound state. This region, however, may already be occupied by nucleons 
of the medium, leading to the Pauli blocking contribution in the quasiparticle shift. Im- 
provements in calculating the wave functions of the isolated nuclei are possible using more 
sophisticated potentials such as BONN, PARIS [2] or AV18/UIX [3] and applying advanced 
methods for the solution of the few-body problem such as the Faddeev-Yakubovski approach 
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19] or the Green's function Monte Carlo method [4j. For tritium, a comparison with 



calculations of the wave function by Wiringa has been performed, and reasonable agreement 
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with the Jastrow ansatz used here has been found. Further improvements of our results for 
the quasiparticle shifts using more advanced approaches to the few-nucleon problem may be 
the subject of future considerations. 



V. APPLICATION: GENERALIZED BETH-UHLENBECK EQUATION AND 
CLUSTER VIRIAL EXPANSION 

The nuclear matter equation of state (EoS) is obtained from Eq. (J3J) after specifying the 
self-energy E(l, z). Taking into account only two-particle correlations in ladder approxima- 



tion <^j, A = 2, we obtain the generalized Beth-Uhlenbeck formula lCj, 111 ] 



n B (T, ii p , fin) = m(T, fjL p , ii n ) + n 2 (T, /i p , /x n ) . (44) 

The single-quasiparticle contribution is n\ = n^ ee + rip ee , where n f * ee (T, /i n , fi p ) = 
2(2n)~ 3 J d 3 P f T (E^ n (P)) describes the free quasi-particle contributions of neutrons and 
protons, see also Eq. (JTj). The two-particle contributions n 2 = n2° und + n| cat contains the 
contribution of deuteron-like quasiparticles (spin factor 3) 

n h r nd (T, » p , Mn ) = 3 / — f d (ET(P; T, /v Mn )) , (45) 

J p > pMott [271)° 

with fd(E) = [e ( - E ~^ ,p ~^ n ^ T — l] -1 , and scattering states of the isospin singlet and triplet 
channel r 2 (degeneration factor j T2 ) 

7^3 J ^ A(A£f (P) + E) sin 2 5%(E, P)^(E, P) . 

(46) 

AE^ E (P) = AEf^P/2) + AEf'{P/2) is the shift of the continuum edge (self-energies 
at momentum P/2), and S^(E,P) denotes the in-medium two-nucleon scattering phase 
shift in the channel t 2 with relative energy E and center of mass momentum P. The 
quasideuteron binding energy E^ U (P;T, fj, p , /^ n ) depends on temperature and nucleon den- 
sities of the medium, expressions are given in Sec. IHIB 21 The shift of the quasiparticle 
binding energy increases with nuclear matter density so that the bound state may merge 
with the continuum of scattering states (Mott effect). P^ ott (T, fi p , /x n ) denotes the momen- 
tum P where, at given temperature and nucleon density, the binding energy of the deuteron 
bound state vanishes. Above the Mott density defined by P^ ott = 0, the integral over P in 
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the bound state contribution n2° und fH5|) is restricted to the region where bound states can 
exist. 

The generalized Beth-Uhlenbeck formula (1441) can be considered as a virial expansion of 



the EoS and was first applied to nuclear matter in Ref. [11] . Results using the PEST4 
interaction potential are shown in Ref. 10j. Due to the inclusion of medium effects such 
as Pauli blocking and self-energy, a smooth transition has been obtained from the low 
density limit, describing nuclear statistical equilibrium between nucleons and deuterons, 
to high densities, where the nucleon quasiparticle picture can be used. Neglecting the 
medium effects, the ordinary Beth-Uhlenbeck formula is recovered, where the single nucleon 
contribution contains the energy E(l) = h 2 P±/2m instead of the quasiparticle energy of free 
nucleons. In the correlated density n 2 , the free deuteron binding energy and scattering phase 
shifts enter the ordinary Beth-Uhlenbeck formula. Here, the term sin 2 5(E) preventing double 
counting if quasiparticle energies are taken in n\, is absent, see lOj. Then, the evaluation of 
the second virial coefficient can be traced back to on-shell properties such as binding energy 
and scattering phase shifts. This has been performed recently by Horowitz and Schwenk 23 1 
using directly the empirical data, instead an interaction potential which is constructed to 
reproduce these data. However, the ordinary Beth-Uhlenbeck formula is restricted to only 
low densities. 

As also shown in Ref. 11] , the contribution of higher clusters to the EoS can be included 



after a cluster decomposition of the self-energy S(l, z), Eq. ([H]). Neglecting the contribution 
of the scattering states, a generalized form of the ordinary NSE is obtained, as given by Eqs. 
( JTOl) . The quasiparticle shift of the nuclei can be calculated in cluster-mean field approxi- 
mation, see Eq. (jl~5l) . Results for light clusters, A < 4, have been given in the previous Sec. 
IIVI Due to the medium dependence of the quasiparticle energies, the contribution of A > 1 
to the EoS (fTUI) fades with increasing density so that a smooth transition from low densities 
to saturation density is obtained. 

The light nuclei have also been considered by Horowitz and Schwenk 23] where a cluster- 
virial expansion was discussed. The second virial coefficient of the virial expansion is gen- 
eralized taking into account also cluster-cluster scattering phase shifts. To obtain rigorous 
results in a given order of density avoiding double counting, and to include medium ef- 
fects, one has to pass over to a systematic quantum statistical approach. However, the 
cluster-virial expansion may be of value when special clusters contribute dominantly to the 
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density. 

In principle, the approach given here can be extended to arbitrary clusters A > 4. Es- 
timates of the quasiparticle shift due to the interaction of the heavier cluster with the sur- 
rounding nuclear matter have been given in Ref. j^j]. For large proton number Z, Coulomb 
effects become important. Calculating the composition of warm dense nuclear matter, the 
mass fraction of nucleons bound in heavy clusters is increasing with decreasing temperature 
and increasing density. This limits the region in the phase space where the EoS is determined 
only by clusters with A < 4. 

The perturbative treatment of the cluster quasiparticle energy shift is justified in the first 
order of density. Higher orders can be determined when comparing with the full solution of 
the medium modified Schrodinger equation for the A-nucleon cluster ( 1151). A term which is 



of second order with respect to the densities has been considered in Ref. [25j. Another effect 
is the formation of quantum condensates due to the Bose distribution function occurring 
in the bound state contribution f j45l) . but also in the scattering contribution fj46|) to the 
density. Degeneracy effects are included in the quantum statistical approach used in this 
work. However, the evaluation of the EoS including quantum condensates such as pairing 
or quartetting needs further consideration. 

We shortly comment the question of thermal stability in connection with the EoS (flOl) . 
It is well-known that symmetric matter at zero temperature becomes instable against phase 
separation below saturation density. The conditions of thermodynamic stability are directly 
related to the behavior of the chemical potential as a function of the densities and will not be 
detailed here. Consequently, the EoS (fTUl) describes nuclear matter in thermal equilibrium 
only outside the instability region. Parameter values belonging to metastable states or 
unstable states may occur in nonequilibrium situations or in inhomogeneous systems when 
a local density approach is considered. They have also to be considered when the Maxwell 
construction is performed to determine the region of instability in the phase diagram. 

The liquid-gas like phase transition in nuclear matter has been considered in the EoS of 
Lattimer and Swesty 2l| and in the EoS of Shen et al. 22|. Including Coulomb interaction, 
an optimum spherical density profile in a Wigner-Seitz cell approach was determined which 
can be interpreted as representative for heavy nuclei, but also for droplet formation. The a 
particle as representative of light clusters has been included. To mimic the density effects for 
the a particle, both approaches used the concept of the excluded volume which cannot be 
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rigorously derived in a quantum statistical approach. In contrast, Pauli blocking considered 
here defines a microscopic process for the medium effects of light clusters. 



VI. CONCLUSIONS 

In the low density limit, the nuclear statistical equilibrium (NSE) with binding energies 
of the isolated nuclei is obtained from the quantum statistical approach to nuclear matter. 
We consider only light clusters, A < 4, but heavier cluster may become more important as 
temperature goes down and density increases. Thus, focusing on only light elements in the 
EoS restricts the temperature and density parameter to the region where the mass fraction 
of heavy elements is small, but our approach may also be extended to heavier nuclei, see 



Ref. 



Deviation from the NSE are due to medium effects, which become relevant once the 
baryon density exceeds 10 -4 fm -3 . A-nucleon correlation effects are described by the A- 
nucleon spectral function which defines the A-cluster quasiparticles. The dependence of 
the cluster quasiparticle energy on temperature and nucleon densities is approximated by 
analytical expressions, Eq. fl33|) for A = 2 and Eq. (j43p for A = 3,4, combined with the 
momentum dependence according to Eq. (|40l) . Compared with more accurate numerical 
calculations, deviations are of the order of 5 %. Analytical expressions for the cluster 
quasiparticle shifts are convenient for calculating the thermodynamic properties of nuclear 
matter in a large parameter range. The evaluation of the cluster-quasiparticle shifts is 
further improved considering more sophisticated potentials and wave functions as obtained, 
e.g., in Green's function Monte Carlo approaches. The values given in the present work are 
approximate estimations, similar as the Skyrme or RMF approaches for the single-nucleon 
quasiparticle case. 

With the shift of the quasiparticle energies, properties such as the EoS ( TTOT) can be deter- 
mined in the subsaturation region. It is possible to interpolate between the low-density limit 
where the NSE is valid, and the saturation density where the single nucleon quasiparticle 
picture can be applied. 

To go beyond the quasiparticle picture, the full A-nucleon spectral function Sa should 
be explored. Instead of 5-like quasiparticle structures, Sa accounts for weakly bound states 
as well as scattering phase shifts including resonances consistently. The full solution of 
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the cluster-mean field approximation would be an important step in this direction. The 
construction of a nuclear matter EoS remains a challenging topic not only in the high-density 
region, see Q], but also in the low-density region where the concept of nuclear quasiparticles 
given here may be a valuable ingredient. 
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APPENDIX A: T-MATRIX APPROACH TO THE SELF-ENERGY 



A detailed derivation of the expressions for the self-energy and the EoS can be found in 
the literature, see, e.g., Refs. Q, Q, Q- 

We give here only some relations, using the 

short notations E x = Ef u (l), E AuP = Ef v (P). 

In Eq. ([8]), the free (A — 1) quasiparticle propagator is 

(o a A- - fi,z 2 (2)-fi,z A (A) 

t^A-llA -> A Z) - —-7 TTT— i rp \ l Ai J 

Z — £ji — ... — &A JA-^Za-iK&Z + ■■■ + &A) 

The T^ matrices are related to the A-particle Green functions 
T A (1. .. A,l'. .. A', z) = V A (l...A,l'... A') 

+V A (1 . . . A, 1" . . . A")G A {1" . . . A", I'" . . . A'", z)V A (l'" . . . A'", I' . . . A') (A2) 

with the potential V A (1 . . . A, 1' . . . A') = J2i<j ^(u> i'f) Ylk^ij ^k,k', an d subtraction of dou- 
ble counting diagrams when inserting the T matrices into the self-energy. The A-particle 
propagator obeys a BSE ([13]) and is solved by the bilinear expansion Q). 

The evaluation of the equation of state in the low-density limit is straightforward. With 

T A {l...A,l'...A',z) = (A3) 
(s - Ei - ... - E A )jj AuP (l . . . A)r AvP {l' ■ ■ ■ A')(E AuP - E v - ... - E A i) 

z- E AvP 

U.P 



24 



we can perform the Q\ summation in Eq. (jSJ). We obtain the result 

y-^ 1 (Q A — Ei — ... — E A )(E AuP — Ei — ... — E A ) _ 

Q\ — z v — E 2 — ... — E A O a — Eaup 



fA-l(E2 + ---+EA) — = — — — ^ /a(Eaup) 



Z v + E2 + ■■■ + Ea — EauP Z u + E2 + ... + Ea — EavP 

Taking ImS(l,z) and integrating the 5 function arising from the pole in the denominator, 
the leading term in density is ji{E Av p - E 2 - ... - E A )f A - 1 {E 2 + ... + E A ) = fA{E A up)- 
Neglecting the contribution of the scattering states, we obtain the generalized form (flQ|) of 
the NSE. 



APPENDIX B: THE CLUSTER-MEAN FIELD APPROXIMATION 



The cluster-mean field approximation [15j is inspired by the chemical picture where bound 
states are considered as new species, to be treated on the same level as free particles. We 
consider the propagation of an A-particle cluster ({A, u, P}) in a correlated medium. The 
corresponding A-particle cluster self-energy is treated to first order in the interaction with 
the single particles as well as with the 5-particle cluster states ({B, u, P}) in the medium. 
The S-clusters in the surrounding medium are distributed according to Eq. ([3]). Full anti- 
symmetrization between both clusters A and B has to be performed, in analogy to the Fock 
term in the single-nucleon case. 

For the A-particle problem, the effective wave equation reads 

A 

[E(l) + ... E(A) - E A up]^Aup(1 ...A)+J2J2 • • • A 1' • • ■ A')^aup(1' ■■■A') 

V.. .A' i<j 

+ E ...A,l>... A')V^p(1' ...A') = , (Bl) 

V...A 1 

with V4{1... A, T... A') = V (12, l'2')5 33 , ...5 aa>- The effective potential 



VJnatter(l . . .A, 1' . . . A') describes the influence of the nuclear medium on the cluster 



rA,mf 
matter 

bound states and has the form 



V££L(1 ...A, T...A') = E A (*)<W • • • <W + E'^ 1 ■ ■ • A 1' ■ ■ • A') . (B2) 



t,3 
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The self-energy like contribution 
A(l) = £v(12,12)ex/(2) 



B 



-EEEE fB,z{ E BvP) £ c 1 • • • ^ i' • • ■ 5')^p(i • • • B)r mP {i' ...&) 

B=2 pp 2...B V...B' i<j 

contains the Hartree-Fock quasiparticle shift. The interaction-like contribution 

AV£(1 . . . A, 1' ... A') = -6w . . . Saa, { I \f(l) + V(12, 1'2') + (B3) 



2 r 

+ f B ,z{ E BoP) £ • • • 1'2 • • ■ ^)fe(22 • • • B)^ BpP (2'2' ...B>) 



B 



B=2 pp 2...B2'...B' 3 

accounts for the Pauli blocking of the potential. The quantity 

oo 

/(I) =A(1) + EEE /b,z(^)I^p(1 • • • ^)| 2 (B4) 

B=2 pp 2...P 

describes the effective occupation in momentum space. The bound states contribute ac- 
cording to their wave function and probability distribution. Note that within the mean-field 
approximation, the effective potential V^tter remains energy independent, i.e. instantaneous. 
Besides the Hartree-Fock term and the Pauli blocking term, determined by the effective oc- 
cupation /(l), the additional terms in AV^ and A(l) account for antisymmetrization. 

APPENDIX C: WAVE FUNCTIONS AND SHIFTS 
1. Quasinucleon 



The nucleon quasiparticle energy shift (1221) contains the nucleon quasiparticle shift AE^ E 
and the effective nucleon mass m*. We expand with respect to the baryon density, 

AE T {P;T,n B ,a) = 5E T (P;T,a)n B + 0(n 2 B ), 

-(T,n B ,a) = l + 6m?\T,a)n B + 0(n 2 B ). (CI) 



I parametrization by Vautherin and Brink [5|, I2JJ which represents an analytical expression 
for the quasiparticle shift of the single nucleon states, 

A£f (0) = | (1 - *„) »n + | (2 + x ) n p+ 1 (2n n n p + r® + (| + ^f) r n + ( | + |) r p , 

(C2) 
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and the effective mass 



m r 



1 + 



*i,3M pi,*2. 

n„ + — n, 

4 4 / " V 2 2 



ft 2 



(C3) 



with the parameter values t = -1057.3 MeV fm 3 , t x = 235.9 MeV fm 5 , t 2 = -100 MeV 
fm 5 , ^3 = 14463.5 MeV fm 6 , xq = 0.56. r n denotes the neutron kinetic energy per particle. 
Expressions for the protons are found by interchanging n and p. For symmetric matter we 
find the shift 5E T (T,0) = — | 1057.3 MeV fm 3 . The temperature dependence is only weak. 

Recently, the EoS of asymmetric nuclear matter has been investigated in the low-density 
region below the nuclear saturation density [ijj]. Microscopic calculations based on the 
Dirac-Brueckner Hartree Fock approach with realistic nucleon- nucleon potentials are used to 
adjust a low-density energy functional. This functional is constructed on a density expansion 
of the relativistic mean-field theory. An improved version of the relativistic mean field 
approach to reproduce properties of nuclei in wide range of A and Z was worked out by Typel 
6J. The relativistic quasiparticle energy ( !23l) reads in the non-relativistic limit AEf E (0) = 
—S(n n ,n p ,T) + V T (n n ,n p ,T) and m*/m T = 1 — S(n n ,n p ,T)/ (m T c 2 ). For the low-density 
limit, expansions of the scalar and vector potentials with respect to the neutron and proton 
densities can be given. 



2. Quasideuteron 

In the deuteron case, we introduce Jacobi coordinates such that q± = {j>2 ~Pi)/2 = 
Pieh 0.2 = Pi + P2 = P or px = —qi + q 2 /2, p 2 = q\ + ^2/2. The kinetic energy is KE = 
T^{2q\ + the interaction is parametrized by the separable Gaussian (TT7I) 

2 / 2 

V(qi, q 2 , q[, q' 2 ) = Ae^e"^^ . (C4) 

The deuteron wave function (1211) results as ip 2 (gi) oc e- 2( tl a2 /{q 2 Jb 2 + 1). The value of the 
rms radius follows as 

rms 2 = -([(n - Rf + (r 2 - R) 2 ]) = ~ J d" qi 

As in the case of single nucleons, the quasiparticle shifts can be expanded as power series 
of the densities, 

AE d (P; T, n B , a) = 5E d (P; T, a)n B + 0(n%) (C6) 
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(C5) 



The first order term SE d consists of the self-energy contribution and the Pauli-blocking con- 
tribution, see Sec. IHIB 2\ in particular ^^ lgldshlft = 5E P (0; T, a)+5E n (0] T, a). Furthermore 
we have 

777* 

— -{T, n B , a) = l + 5m d (T, a)n B + 0(n%) . (C7) 

Values for 5E^ uli (0; T) are given in Tab. HV1 For T=10 MeV we have 8m d (10, 0)=21.3 
fm 3 , whereas for T=4 MeV the value 8m d (A, 0)=87.1 fm 3 results. Due to the Pauli blocking 
both quantities are strongly temperature dependent. At zero temperature and low densities, 
we find for the Gaussian interaction 

SE d (P,0)---A^(P/2)e J d 3 qilMqiW ■ ( C8 ) 

3. Quasitriton/helion 

Next we consider A — 3 (t, h). Jacobi coordinates are q\ = \{p2 — Pi), 92 = §( — §pi — 
\P2 + P3), 93 =Pi+P2+P3 or pi = -qi - \q 2 + §93, P2 = qi - |?2 + §g3, P3 = Q2 + |<?3- The 
kinetic energy is KE = ^(2qf + \q\ + 

We start from a Gaussian pair interaction (fTTI) which gives in Jacobian coordinates the 
three-nucleon interaction 

2 ; 2 



^"(gi, 92, 93, 9l, 9 2 > 9 3 ) = j<W e ^ e 

(qi+|g2) 2 (n-^n+^ij) 2 (gi-jm) 2 (gi+^272^) 2 ^ 

+5, ,_ ffi e e ^ +5, j e ^ e ^ k (C9) 

The Gaussian variational ansatz (fl8"|) reads after introducing Jacobians (indices in B, a, b 
are omitted) 

^ 3 GauSS (<7i, 92, 93) oc e-S e -|#5 93i p . (CIO) 



The Jastrow variational ansatz (121]) motivated by the solution of the two-particle problem, 
reads after introduction of the reduced Jacobian coordinates Xi = qi/b and choosing the 
coordinates as x± = £i{(l — z 2 ) 1 / 2 , 0, z}, x 2 = X2{0, 0, 1} 

_ 3 i> 2 T 2_ 9 6 2 9 

Jastrow e 2 a 7 8 ^ 

^3 a ™(x 1; X 2 ,z) OC (x 2 + 1)( l x? + A X 2 + l XiX2Z + 1)( 1 X 2 + _ |^ + fj ( CU ) 

The kinetic energy follows as 

KE 3 = f dx2x\ [ dx\ x\ ( x\ + —x 2 , J / dz<p\(xi, x 2 , 2) (C12) 

™ A^ 3 J J V 4 / 7-1 
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with the norm 



N 3 = I dx 2 x\ I dx\%\ I dz(fil(xi, X2, z) 



For the potential energy we obtain 

b 3 



PE« = 3A, 



dx 2 X2 



_ 2 b A 

dx\x\e Xl ~? j dz(fs(x\, x 2 , z) 



The nucleonic point rms radius follows as 

rmso = — — / dx 2 xl / dx\ x] / dz 
3 b 2 N 3 J 2 J V-i 



-1 2 



(C13) 



(C14) 



1 ( d(p 3 \ 2 ^ 2 ( <9y? 3 x 1 



6 V dxi J 9 V dx 



(C15) 



in particular rms^ = 2/B 2 for the Gaussian ansatz (jl8p . 

At finite temperature, the Pauli blocking contribution to the quasiparticle shift is given 

by 



5£ 3 Pauli ' J (P) 



3A,6 3 f2nh 



]ir 2 N 3 V TflT 



2\ 3/2 



(C16) 



^2 ^2 

x / dxix\ I dzi<p 3 (xi, x 2 , z 1 )e~ XlT i 7 e~^T^ Xl+ ^ X2 ~^ X3 ' )2 I dx^x\ I dz 5 ip 3 (x 5 , x 2 , z 5 )e~ XsT > 7 



,2 if 



At zero temperature where p\ ~ (xi = x 3 /3 — x 2 /2; x 3 = 0), the shift is given by 
3A 3 



dx 2 x 2 - 



2N 3 



_ X .2(|4 + J>L) 



(^ + l)e + l) 2 



/ dz 5 ip 3 (x 5 ,x 2 ,z 5 )e Xs ^ . (C17) 



4. The a-quasiparticle 



To solve the four-nucleon Schrodinger equation in the zero-density limit, we separate 
the center-of-mass motion from the internal motion introducing Jacobian coordinates, q x = 

-|pi + §P2, Q2 = -5P1 - 5P2 + §P3, 93 = -jPl - \P2 ~ jP3 + |P4, ?4 = Pi + P2 + P3 + Pi- 

The inverse transformation is p x = \q A - \q 3 - \q 2 - q x , p 2 = \q^ - \q 3 - \q 2 + qi,p 3 = 

IQa ~ §<?3 + ?2, Pa = \q± + Q3- 

The Schrodinger equation separates with <y9 Qj p(l, 2, 3, 4) = ^((fi, g*2, <?3)<5p <? 4 - The kinetic 



29 



energy is KE4 = ^(2g 2 + |gf + |gf + \P 2 )- The potential energy follows as 



{2 / 2 

(91 + |g2) 2 (gl-^g2+292) 2 (gl-|g2) 2 (91 + ?92 — ggg > 2 



(-92 + f 93) (-3q 2 -2q 3 +6^) 2 



+ ^l, 9 ' 1 ^, g 2 + |fe- g 3) e 6 

(gl+jgg+j 93) 2 (6 gl +3 g2 -4q3 + 12^) 2 

5 fe-9 3 ) 92-92-3(93-93) 

(-gl + ^g2 + §93) 2 (-6 9 i+3.; 2 -493 + 12 9 3) 2 
+ 5 9i,9 1 + |(93-9^)^,92-i( g 3-9^) e ™ e ^ I ^ ( C1§ ) 

To solve the internal motion we use a variational ansatz for the wave function. The 
Gaussian variational ansatz (jT8l) reads after introducing Jacobians 

2q 2 3q% 4qf 

^ auss (gi, g 2 , g 3 , ? 4 ) oc e ^5 g4>P . (C19) 

The Jastrow variational ansatz (1211 



7^(P2-Pl) 2 -i^(P4-p 3 ) 



^t w (pi, a, Pi) e (P2 4 : P1)2 , - • • • e (P4 4 : P3 ) 2 , , w**^ • < ( ,2 ° ; 



4b 2 ^ x 46 2 

reads after introduction of the reduced Jacobian momenta Xi = (ji/b, with x% 
■Ti{a/1 - z\ cos(0i), a/1 - i? sin(0i), ^j; f 2 = z 2 {0, 0, 1}; x 3 = £ 3 {a/1 - if, 0, z 3 }, 

g _ ^( 2:E l + | x2 + | x 3) 



+ l)[(\x\ + + l) 2 - ^x\xlzl]{\xl + fx 2 - \x 2 x 3 z 3 + 1) 

1 



[(\ X l + l^ X 2 + 1^3 + |^2^3^3 + I) 2 - {\X\X*Z X + |xiX 3 {2i2 3 + a/ 1 - z\ a/ 1 - if COS0i}) 2 ] 

(C21) 



We evaluate the norm as 



2 Z" 1 d2! f™ 2 f™ 2 f 1 dz 3 2 



iV 4 = I dx x xi J -j- J J dx 2 xi J dx 3 xi J . (C22) 



The kinetic energy of the internal motion is calculated from 

KE4 = Si dxixl S-^L dx2xl l dX3xl L 
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The potential energy is (all six terms give the same contribution; we take the first one as 
representative) 

2 



PE 4 



6A6 3 
2tt 2 N a J 



1 dz 3 ( f°° , 2 f 1 dz x r2a 



dx 2 x 2 / dx 3 x 3 / — - I / dx 1 x 1 / — - 



i 2 



i 2 



2tt 



(C24) 



The integrals over <\>\ can be performed analytically. 
The nucleonic point rms radius follows as 



rms. 



/ dx 1 x 1 / — / — / dx 2 x 2 / dx 3 x 3 / — 

t/Q J — 1 *-* »/ « — 1 



4 Jo 
1 / <9y9 4 



dxix\ J 

2 ! / o \ 2 



1 / 8(f 4 



8 \dxi J 6 \(9x 2 



3 / d(f4 



16 V <9x 



(C25) 



in particular rms 2 = 9/4B 2 for the Gaussian ansatz (fT8|) . 

Calculating the rms radius for the Jastrow ansatz, the three terms give the same contri- 
bution so that we take the last one as representative 



rms. 



x < Xn 



16 b 2 N d 



r a 2 r ^ r d ^ r a 2 r * 2 f 1 

i dx\X l I / / dx 2 x 2 / dx 3 x 3 I 

o J-i 2 Jo 27r J J • J_ x 



1 d,Z3 2 

2"^ 4 fhfi - fi) 2 



£ ( £2 £2\ / £l £l\ £ £ 

^2-/11/2 _ J3J + "q"v/2 ~~ h) + -ghh 



+ x 2 2 



"3(^2 — /I) + 3A/2 



- 2 




+ x 2 


— T^JiJa 



3^/i(/ 2 2 -/3 2 ) + ?(/ 2 2 -/3 2 ) + ^/i/2 



9 



+2XiX 3 (ZiZ 3 + ^/l - ^ 2 ^/l - 2§C0S(^ 

+2a;ia;2^i 



9 



"g(/ 2 2 _ /I) + 3/1/2 



THphifi ~ fi) + If (/a ~~ /I) + "^"/i/ 2 



16 ff 



"3I/2 ~ /a) + 3/1/2 



16 f f 

-g-71/3 



(C26) 



where we used the abbreviations /x = + |x 2 — \x 2 x 3 z 3 + 1), / 2 = (|x 2 + ^x 2 + + 
|x 2 x 3 2;3 + l),/ 3 = (\x x x 2 zi + \x\x 2 3 {ziz 3 + a/1 - z[^Jl - 2;|cos0i}). 

At finite densities, we take into account self-energies and Pauli blocking due to the 
medium. We obtain 12 terms which differ only by the isospin (p, n) dependence so that 

A £Pauii (p) = _ 6 < P (l,2,3,4)[/ p (l) + / n (l)]K(l,2;l',2')^,p(l / ,2',3,4). (C27) 

1234,1'2' 

First we consider P = 0. At T = 0, the Fermi distribution function is replaced by the step 
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function. In the low-density limit where p\ — or x± — — \x 2 — \x*z, we have 

^Pauli, J(() . Q) = __6A_ ^ ^3 j~ ^ j~ ^ 



X 



{xl + |x§ + \x 2 x 3 z 3 + l)[(f a?| + ±x§ + \x 2 x 3 z 3 + l) 2 - %\x\ + |x 2 x 3 z 3 ) 2 ] 

1 

X 

{A + f ^3 - |^2X3^3 + + ¥ X l + ^ 2 X 3 Z 3 + l) 2 - {\x\ + §X 2 + \x 2 X 3 Z 3 f\ 

x <J dx 5 x 5 J — J — e 1 2 5 yi(x b ,x 2 ,x 3 ,z 5 ,z 3 ,(t) 5 )> • (C28) 

The integral over 05 can be performed analytically. 

For arbitrary temperatures we find in the low-density limit, where we have f T (p) = 
|n r (27r/i 2 /mT) 3/ ' 2 exp[— 2^p 2 ], the expression 

snT) = (^f j- dxsxl jy X2xl (C29) 

- 1 d Zl r j 2 f 2?v #5 



X 



2tt 



x^ 4 (xi,x 2 ,x 3 , 2:1,2:3,01)} I y Gfe 5 x 5 y — y — e t 2 5 </?4 (a*, x 2 , x 3 , z 5 , z 3 , 5 

For finite P we have to introduce in the Fermi distribution pi/b = ^P/b— |x 3 — |rr 2 — X\. 
Expansion for small P gives the Pauli blocking contribution to the effective mass. 
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